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&BS  TRA  C' 


A  global,  six-level,  primitive  equation  model  was  numer¬ 
ically  integrated  to  examine  the  interaction  of  a  aid-lati¬ 
tude  baroclinic  wave  with  several  different  mountain-valley 
configurations.  The  purpose  was  to  determine  the  contribu¬ 
tion  of  a  well-developed  baroclinic  wave  to  the  initiation 
of  a  monsoon  surge.  It  was  found  that  these  waves  can  ini¬ 
tiate  a  surge.  Such  surges  are  weak  and  limited  in  extent 
but  they  do  exhibit  the  characteristics  of  observed  monsoon 
surges.  The  latitudinal  location  of  the  east-west  mountain 
range  has  a  large  impact  on  the  strength  of  the  disturbance 
and  thus  on  the  initiation  of  the  surge.  However,  the  size 
of  the  valley  is  not  an  important  consideration.  The  surge 
is  strongest  at  the  lowest  levels  and  dissipates  rapidly 
with  height.  Overall  results  indicate  surge  forcing  by  the 
baroclinic  wave  is  possible  but  other  forcing  mechanisms  are 
necessary  in  order  to  simulate  the  stronger,  more  extensive 
surges  observed  during  the  winter  monsoon. 
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INTRODUCTION 


In  recent  years  meteorologists  have  given  considerable 
attention  to  the  Asian  winter  monsoon.  Described  by  Chang 
and  Lau  (1980)  as  "one  of  the  most  energetic  heat  engines  in 
driving  the  earth*s  atmosphere",  the  winter  monsoon  is  the 
most  distinctive  component  of  the  Northern  Hemispaere  winter 
circulation.  Its  fundamental,  large  scale  structure  has 
long  been  observed,  however,  much  about  t  h  .or.soon  remains 
unknown.  The  role  of  the  monsoon  in  the  :  traction  of  the 
mid-latitude  and  tropical  regions,  the  v  '  ility  of  the 
monsoon  and  the  forcing  cf  small  scale  monsoonal  variations 
are  all  areas  of  active  research. 

A  mcnscon  may  be  defined  as  a  three  dimensional,  plane¬ 
tary  scale  wind  regime  that  reverses  direction  seasonally. 
In  the  case  of  the  monsoon,  the  establishment  of  the  summer 
cyclone  and  winter  anticyclone  is  governed  by  land-sea  dif¬ 
ferences  and  the  varying  position  of  the  sun.  Essentially, 
it  is  the  seasonal  replacement  of  one  persistant  circulation 
system  by  an  equally  persistant  but  opposite  system.  The 
development  of  these  circulations  closely  links  the  monsoon, 
pa rt icular ily  that  of  east  Asia,  to  regional  topography. 
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The  long,  clear  nights  of  the  Siberian  winter  allow 
strong  radiaticnal  cooling  and  the  subsequent  formation  cf  a 
shallow,  cold  layer  of  surface  air.  Persistent  cold  air 
advection  throughout  the  troposphere  further  contributes  to 
the  formation  of  a  surface  high.  The  winter  monsoon,  which 
dominates  the  entire  Asian  continent,  is  characterized  by 
this  intense,  cold  polar  anticyclone.  The  resulting  north- 
northeasterly  surface  flow  extends  over  all  of  eastern  Asia 
and  southward  across  the  equator.  The  Himalayan  Mountains 
restrict  the  southward  movement  of  the  cold  lower  tropos¬ 
pheric  air  and  provide  an  effective  barrier  to  the  mixing  of 
the  cold  continental  airmass  with  the  warm  tropical  airmass 
of  southern  Asia.  As  the  anticyclone  becomes  established 
over  Siberia,  a  surface  cyclone  develops  over  the  "maritime 
continent”  region  of  Indonesia  and  Malaysia  as  a  major  plan¬ 
etary  scale  convective  area  is  shifted  eastward  from  its 
summer  position  near  India  (Hamage,  1971)  . 

An  analysis  cf  the  monthly  mean  velocity  potential  field 
at  200  mb  for  December  1  9  74  by  Char.g  and  Lau  (1  980)  shows  a 
broad  region  cf  divergence  over  the  maritime  continent. 
This  is  indicative  of  deep  rising  actions  and  large  scale 
overturning  of  the  atmosphere.  Their  analysis  also  shows  a 
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region  of  maximum  convergence  over  nortaern  China  in  the 
vicinity  of  the  Siberian  anticyclone.  (Note  that  for  large 
scale  flow  the  divergence  pattern  is  approximately  given  by 
the  velocity  potential.)  Upper-level  convergence  is  also 
evident  in  the  equatorial  west  Pacific  Ocean  and  over  the 
east  African  coast.  Although  much  weaker,  these  latter  two 
areas  indicate  the  existence  of  a  talker  circulation. 
Ascending  over  the  naritima  continent,  one  branch  flows 
westward  across  the  Indian  Ocean  and  Arabian  Sea,  while 
another  branch  flows  eastward  to  descend  over  the  western 
Pacific  Ocean. 

Thus  the  winter  tropospheric  circulation  in  this  region 
is  dominated  by  a  strong  local  Hadley  cell;  a  thermally 
direct  circulation  involving  descent  in  the  cold  air  over 
the  northern  anticyclone,  southward  surface  flow,  ascent  of 
warm  air  over  the  maritime  continent  and  northward  return 
flow  at  200  mb. 

The  large-scale  components  of  the  monsoon  are  consistent 
and  identifiable  from  year  to  year,  but  they  are  not  static. 
The  primary  features  not  only  vary  in  strength  and  location 
on  an  annual  time  scale  but  there  are  small  time  scale  vari¬ 
ations  in  the  velocity,  temperature  and  pressure  fields  as 
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The  monsoon  surge  is  initiated  'as  the  upper-level  trough 


moves  eastward  and  the  pressure  gradient  across  the  east 
China  coast  tightens.  The  normal  geostrophic  balance  of  the 
anticyclonic  flew  is  broken  since  the  pressure  gradient 
force  is  now  greater  than  the  Coriolis  force  and  a  rapid 
cross-iso  baric  flow  from  high  to  low  pressure  begins.  This 
accelerating  flow  towards  lower  pressures  at  lower  tropos¬ 
pheric  levels  is  the  monsoon  surge.  Following  the  surge,  an 
increase  in  the  northeasterly  flow  can  be  traced  to  the 
equatorial  regions  within  12-24  h.  Within  one  to  two  days 
convection  in  the  equatorial  region  intensifies  and  sustains 
the  already  enhanced  Hadley  circulation.  The  upper-level 
outflow  from  the  equatorial  South  China  Sea  accelerates  not 
only  northward,  but  within  one  day  of  the  surge  accelerates 
east  and  west  along  the  equator,  thus  strengthening  the 
Walker  circulation.  The  surge  ends  as  the  mid-latitude 
trough/ridge  pattern  moves  sufficiently  eastward  so  that 
subsidence  no  longer  occurs  over  the  anticyclone. 

Lim  ar.d  Chang  (1981)  and  Chang,  Chen  and  Hillard  (  1  983) 
suggest  that  the  monsoon  surge  is  char acterized  by  two 
stages  separated  by  an  interval  of  a  few  hours  to  one  day, 
depending  on  the  location  of  the  observing  station.  The 
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first  stage  is  the  pressure  surge,  or  Leading  edge*  of  the 
downward  accelerating  air,  with  a  propagation  speed  charac¬ 
teristic  cf  gravity  waves.  The  second  stage,  which  moves  at 
advective  time  scales,  is  due  to  a  frontal  passage  and  is 
characterized  by  a  sharp  decrease  in  the  surface  dewpoint. 

Few  studies  have  been  done  to  determine  the  cause  or 
nature  of  the  forcing  that  precipitates  a  monsoon  surge. 
Lim  and  Chang  (1981)  used  the  linearized  shallow-water  equa¬ 
tions  on  an  equatorial  beta-plane  in  an  analytical  study  of 
monscon  response  to  a  transient  mass  source  that  simulated 
intense  anticyclogenesis.  Their  results  showed  that  such 
mid- latitude  forcing  excites  mainly  Rossbv  waves  after  an 
initial  period  of  gravity  wave  type  motions  associated  with 
strong  northerly  winds.  They  suggested  that  pressure  forc¬ 
ing,  using  equatorial  beta-plane  dynamics,  can  explain  the 
large  features  of  synoptic  scale  motions  in  the  eastern  mon¬ 
soon  region.  However,  they  do  not  consider  other  effects 
such  as  orography  or  a  variable  mean  zonal  wind.  Nor  does 
their  method  account  for  the  shallow  vertical  structure  of 
monsoon  surges  or  upper  tropospheric  response  to  the  surge. 

The  purpose  of  this  study  is  to  examine  the  interaction 
of  topography,  specifically  a  high  east-west  mountain  range. 
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with  a  mid-latitude  baroclinic  wave.  The  purpose  will  be  to 
determine  the  contribution  of  the  baroclinic  wave  *.c  the 
initiation  of  monsoon  surges,  k  global,  six-level  primitive 
equation  model  will  be  used  to  simulate  the  meteorological 
conditions  prior  and  subsequent  to  surge  initiation.  Ini¬ 
tial  conditions  will  include  a  strong  baroclinic  jet  north 
of  the  mountain  range  and  a  small  disturbance  in  the  mean 
current;  the  disturbance  will  grow  through  baroclinic  insta¬ 
bility.  Several  mo untain -valley  configurations  will  be  used 
to  study  the  surge  and  numerical  solutions  will  be  compared 
with  observed  monsoonal  surges. 
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II.  DESCRIPTION  OF  MODEL 

The  numerical  model  used  in  these  experiments  is  a  ver¬ 
sion  of  the  UCLA  potential  enstrophy  conserving,  general 
circulation  model  for  an  adiabatic,  dry  and  frictionless 
atmosphere.  The  governing  equations  for  this  model  are  the 
primitive  equations:  the  horizontal  momentum  equation, 

thermodynamic  energy  equation,  surface  pressure  tendency 
(continuity)  equation,  and  the  hydrostatic  equation.  These 
are  given  in  Appendix  A.  The  design  and  dynamic  processes 
of  the  model  are  fully  discussed  by  Arakawa  and  Lamb  (1977, 
1981) .  1  summary  of  the  version  used  for  this  experiment  is 

presented  here. 

A.  THE  HORIZONTAL  GRID 

The  prognostic  variables,  »ind(U,Vl,  temperature  (T)  , 
and  terrain  pressure (n) ,  are  staggered  horizontally  accord¬ 
ing  to  Arakava's  Scheme  *C*  (Fig.  1).  This  grid  arrange¬ 
ment,  although  using  more  computer  time  than  an  unstaggered 
grid,  gives  a  better  structure  for  the  shorter  waves  that 
arise  due  to  geostrophic  adjustment,  nonlinear  interaction 
or  mountain  effects.  Resolution  is  2.5°  long  by  4.0°  lat 
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using  a  36  x  46  spherical  grid.  The  computational  domain  is 
a  90°  sector  of  the  globe  with  cyclic  continuity  imposed  at 
the  east  and  west  boundaries. 

In  spherical  coordinates  the  poles  are  singular  points. 
Since  velocity  components  cannot  be  defined  at  the  poles, 
they  are  taken  as  rr-points.  Due  to  the  convergence  of  the 
meridians  towards  the  poles  a  very  short  time  step  would 
normally  be  required  to  insure  computational  stability.  To 
allow  a  longer  time  step,  longitudinal  averaging  is  done  for 
selected  terms  in  the  equations  that  govern  gravity  wave 
propagation.  To  do  this  the  zonal  pressure  gradient  and 
zonal  mass  flux  are  expanded  into  a  Fourier  series  around 
each  latitude  circle,  then  the  amplitude  of  each  wave  compo¬ 
nent  is  reduced  by  the  factor  S.  This  technique  smooths  the 
gradients  of  the  meteorological  fields  rather  than  their 
Fourier  expansions.  The  smoothing  operator  S  is  applied 
only  at  higher  latitudes  (where  S<1)  and  is  given  by: 

s  _  AX.  £2S$. _ 

Ap  sin(nAA/2) 

(Symbol  definitions  are  given  in  Appendix  3( 
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THE  VERTICAL  GRID 


The  equations  used  in  this  model  are  written  using  the 
vertical  sigma  coordinate: 

a  ■  P  Pt  where  II  3  p  - 

_____  s  t 

P5  is  the  surface  pressure  and  ?t  is  the  pressure  at  the 
tropopause . 

The  vertical  domain  extends  from  the  earth's  surface  to 
200  mfc.  Variables  are  staggered  vertically  (Fig.  2)  in  six 
layers  equally  spaced  in  sigma.  The  lower  boundary,  <7=1, 
follows  topography,  where  the  surface  pressure  is  defined. 


At 

both 

the  upper 

and  1 

ow 

er  b 

oun 

da 

ries 

ver 

tical 

veloc 

1.  * 

y. 

<y. 

va 

nishe; 

s.  Th 

e  pr 

ognost 

ic 

var 

iab 

le 

s  (u 

,  v,  r 

,) 

ar 

3  comp 

ut 

ed 

at 

mi 

d-lay 

er ,  wh 

ile 

geo  pot 

en 

tial 

•  2 

t 

vert 

ical 

velo 

city. 

a, 

and 

1 3 

rrain 

press 

ure. 

n  are 

com 

put 

ed 

di 

agno 

st: 

Lea 

ily  on 

si 

gma 

la 

vels. 

c. 

FINITE  DI 

FFERENCING 

S 

CHEN 

ES 

The 

time 

diff 

er  e  nci 

ng 

met 

hod 

- 

used 

in 

th: 

Ls 

model 

ar 

a 

the 

f  0 

urth- 

order 

space  di 

ff 

eren 

c  in 

g 

sch 

erne 

a: 

id 

the 

Ha 

ts 

uno 

sc 

heme 

(also 

call 

ed  t  he 

E 

uler 

-ba 

ck 

ward 

sen 

em; 

i)  . 

One 

Ma 

ts 

uno 

19 


i 


time  step  is  performed  followed  by  four  centered  time  steps. 
The  Matsunc  scheme  is  used  to  eliminate  the  computational 
mode  of  the  centered  scheme  and  because  it  selectively  damps 
high  frequency  waves.  The  time  step  used  hare  is  360  s. 
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III.  MODEL  INITIALIZATION 


The  modal  uses  an  analytical  initialization  scheme 
developed  by  Major  J.  L.  Hayes,1  OSAF.  A  small  perturba¬ 
tion  is  superimposed  on  a  horizontally  and  vertically  vary¬ 
ing  mean  current,  then  allowed  to  grow  through  baroclinic 
instability.  The  scheme  was  derived  specifically  for  flow 
over  a  flat  earth;  initialization  over  high  terrain  requires 
an  adjustment  to  the  scheme  and  is  discussed  separately. 

A.  INITIALIZATION  OF  THE  MEAN  FLOW 

A  baroclinic  Eicklay  jet  is  situated  at  U5°N,  where  it 
reaches  a  maximum  velocity  of  46  m/s  at  200  mb.  The  initial 
mean  surface  flow  has  the  same  latitudinal  variation  as  the 
jet;  it  reaches  a  maximum  of  5  m/s  at  43°N  and  decreases  to 
zero  at  the  poles.  The  v -component  of  the  wind  is  initially 
set  to  zero.  The  surface  wind  is  described  by  the  Bickley 
jet  profile; 

U  ($)  *  U  sech2  (y  ) 

S  IB 

v  ($)  =o 

s 


1  Personal  communicat ion ,  present  affilitation  Detachment  1, 
HQ  AWS. 
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The  vertically  varying  tern  T  (p ,p)  is  obtained  from  the 
NACA  standard  atmosphere: 


TCp,4>o) 


=  t(pq,$q)  -  rz 


where :  Z  = 


The  constants  are  defined  in  Appendix  3. 

3.  I NIT I ALIZ ATICN  07 ZR  HIGH  TERRAIN 

The  mountains  in  this  model  are  intended  to  be  a  very 
simple  representation  of  the  Himalayan  Range-  The  mountains 
extend  east  to  west  across  the  grid  section  with  a  maximum 
elevation  of  3000  m.  The  longitudinal  spar,  is  16°  or 
approximately  1600  km. 

The  mountains  were  initialized  at  full  height  without 
any  growth  period.  The  equations  used  to  produce  the  moun¬ 
tains  are  relatively  simple  and  produce  a  knife-edged  range: 


$  =  Cl(<j))'Jf  = 


.25  t  at  j  |  *  2 A4> 

.75  *  at  | i  *  2A$ 

<>f  at  $>  = 


The  immediate  initialization  of 
tne  gecstropnic  adjustment  time. 


the  terrain  doss  lengthen 
However,  the  model  is  run 
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for  72  h  prior  to  creating  the  valley  and  it  is  reasonable 
to  assume  the  longer  adjustment  period  will  have  no  serious 
effect  on  the  attempted  simulation  of  the  monsoon  surge. 
Using  the  mean  flow  initialization  scheme  without  a  mountain 
growth  period  resulted  in  a  large  scale,  north-south  oscil¬ 
lation  in  the  pressure  field.  Thera  was  little  evidence  of 
damping  during  a  95-h  model  run.  To  solve  this  problem  the 
surface  pressures  at  mountain  grid  points  were  adjusted  for 
height  above  mean  sea-level.  Using  the  hydrostatic  equation 
and  the  adiabatic  lapse  rate  for  the  NAUA  standard  atmos¬ 
phere,  the  following  equation  was  developed  for  a  new  (moun¬ 
tain)  surface  pressure  in  terms  of  the  old  (flat  earth) 
surface  pressure: 

/  rz\ g/RT 

P  (new)  =  P  (old)  (1  -  —  j 

s  s  y  tq  j 

The  new  surface  pressure  is  used  to  re-calculate  wind  and 
temperature  values  at  each  mountain  grid  point.  No  oscilla¬ 
tion  was  evident  when  the  adjusted  values  were  used  in  the 
initialization. 

As  the  monsoon  surge  is  expected  to  occur  almost  immedi¬ 
ately  following  the  creation  of  the  valley,  the  geostrophic 


26 


Thre¬ 


ad  j  ustmant  time  cannot  ba  neglected  at  this  point, 
fora,  the  valley  is  "built1*  ovar  a  12-h  period  to  minimize 
the  tiae  required  for  gaostrophic  adjustment.  The  method 
used  is  based  on  the  scheme  developed  by  Hayes  (previously 
cited)  for  raising  terrain  height  during  a  model  run.  To 
build  the  valley  however,  the  method  is  reversed: 


«  = 


Cj  (*)4f  £  1  -  C2  ( t )  G ( X ) J 


?he  function  C,(r)  is  defined: 


C2(t)  = 


sin2|  ¥0 


Ti<T<Tf 


T4T.  and  t>t 
i  f 


The  function  G(,\)  is  defined: 


GM  = 


X-Xi-2AX 

-2AX 

X?  +2AX-X 
2AX 

0 


Xj<X<X2 
(X-2AX)4 X*Xi 

X2*X*(X2+2AX) 

elsewhere 


\,  and  are  the  west  and  east  edges  of  the  valley  respec¬ 
tively.  Where  no  valley  is  built,  G(X)  =  i.  The  valley 
walls  are  sloped  over  two  grid-lengths  ( 2 AX)  ,  and  the  valley 
floor  is  at  mean  sea-level  (0  m)  . 
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C.  THE  PERTUREATI3N 

The  perturbation  is  a  standard  barotropic  disturbance. 
It  consists  of  a  weak  longitudinal  sine  wave  with  maximum 
latitudinal  amplitude  at  45°N.  The  vertical  amplitude  is 
one-tenth  the  earth's  radius.  The  perturbation  is  described 
by : 

#’  =  fg  (A  sin(nX)  sin2  (2$)) 

The  following  perturbation  quantities  are  derived  using  the 
hydrostatic  equation  and  geostropnic  wind  equation: 

P‘  •  -  H  =  -  H°(a  sin(nA)  sin2  (2+)) 

n.  *  =  -  ifA  sin(nX)  2sin(4<t>)) 

afgdtp  av  / 

_  1  6 <t> *  _  1  fnA  cos(nX)  sin2l2iM) 

~  afo  cos  5  4>  a  cos$  '  / 

T’  =  0 

The  total  fields  are  calculated  by  adding  the  perturbation 
quantities  to  the  lean  current.  The  total  fields  will  then 
be  integrated  ir.  time  to  arrive  at  the  required  forecast 
values. 
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IV.  ANALYSIS  OF  RESULTS 

A.  EXPERIMENTS  152:  CONTROL  RUNS 

Two  control  runs  were  mads  for  comparison  with  subse¬ 
quent  model  runs.  The  first  run  uses  a  flat  earth  while  the 
second  includes  a  3000-m  east-west  mountain  range  centered 
at  3 0°N.  Both  runs  are  initialized  with  the  barcclinically 
unstable  mean  flew  and  the  barotropic  disturbance  discussed 
in  Chapter  3.  Model  times  given  are  based  on  the  time  (h) 
at  which  the  disturbance  is  introduced  into  the  mean  flow. 

The  first  run  follows  the  development  of  the  disturbance 
independent  of  terrain  effects.  Sea-level  pressure  fields 
for  60,  72,  84  and  96  h  are  shown  schematically  in  Figs.  3 

and  4.  Potential  temperature  fields  for  the  same  times  are 
shown  in  Figs.  5  and  6.  Since  a  pie-shaped  section  of  the 
globe  is  represented  in  cartesian  coordinates  there  is  some 
distortion  of  the  fields.  This  is  primarily  evident  in  the 
tight  north-south  pressure  gradients  vaioh,  for  comparable 
wind  velocities,  appear  weaker  in  the  east-west  direction. 

The  disturbance  is  strictly  limited  to  the  mid-latitudes 
with  a  very  fiat  pressure  gradient  in  both  the  polar  and  the 
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tropical  regions.  The  disturbance  moves  eastward  at  a  con¬ 
stant  velocity  of  approximately  8.5°  every  twelve  hours. 
Wave  development  is  realistic  with  a  995  mo  low  and  a  1023 
mb  high  developing  after  78  h. 

Development  of  the  disturbance  in  the  wind  and  tempera¬ 
ture  fields  closely  parallels  that  of  the  pressure  field, 
although  the  temperature  and  potential  temperature  data  dis¬ 
play  far  more  short  wave  noise.  As  no  surface  frictional 
dissipation,  diffusion,  or  smoothing  techniques  are  used  in 
the  model,  some  noise  in  the  results  is  expected.  The  large 
scale  temperature  features  remain  discernaole  throughout  the 
six  day  (144  h)  run  but  the  the  small  scale  features  are 
masked  by  noise  after  126  h,  especially  at  higher  levels. 
(Temperature  and  wind  velocities  are  plotted  on  sigma  sur¬ 
faces;  'higher  levels'  refers  to  increased  height,  or  lower 
sigma  values.) 

After  96  h  particularily  strong  cyclogenesis  begins  to 
occur.  At  114  h  the  cyclones  intensify  rapidly  and  the 
ridge  is  drastically  weakened.  While  diffusion  or  smoothing 
applied  to  the  fields  would  decrease  the  strong  gradients 
caused  by  cyclogenesis,  they  would  also  weaken  the  gradients 
of  pressure  and  temperature  which  could  contribute  to  possi- 
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bis  monsoon  surges.  For  this  reason  neither  method  is  used. 
k  linear  drag  tern  was  applied  to  the  winds  in  the  lowest 
level,  however  it  had  no  discs rnable  effect  on  the  strong 
cyclogenesis. 

The  mountain  range  used  in  second  control  run  is  the 
basis  for  all  subsequent  terrain  configurations.  Fig.  7 
gives  a  three-dimensional  raprssentation  of  the  range. 
(Note  the  vertical  scale  is  greatly  exaggerated. )  The  data 
from  this  run  are  used  to  determine  the  time  and  location  at 
which  the  valley  is  to  be  created.  The  requirements  are: 
(1)  a  high  pressure  center  to  the  northeast  (to  be  consis¬ 
tent  with  the  Siberian  anticyclone)  and  (2)  a  sea-level 
pressure  difference  of  25  mb  between  the  high  and  low  pres¬ 
sure  centers.  The  latter  requirement  ensures  an  adequate 
northwesterly  current  to  initiate  a  mor.soonai  surge. 

Sea-level  pressure  and  potential  temperature  fields  for 
60  ,  72,  84  and  96  h  are  shown  in  Figs.  3  through  11.  East¬ 
ward  movement  of  the  disturbance  is  again  3.5°  per  twelve 
hours.  The  most  noticeable  difference  between  this  experi¬ 
ment  and  the  first  is  a  series  of  weak  highs  and  lows  that 
have  developed  along  the  northern  slope  of  the  mountains. 
These  disturbances  migrate  eastward  at  a  spaed  comparable  to 
the  mid- latitude  wave. 
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Development  of  the  wave  is  slower  than  in  the  previous 


case.  At  78  h  a  999  mb  low  and  a  1019  mb  high  have  devel¬ 
oped,  giving  a  pressure  difference  of  20  mb  compared  to  28 
mb  earlier.  This  represents  an  overall  weakening  of  the 
disturbance  and  is  not  due  to  the  weakening  of  either  the 
high  or  low  centers  individually.  It  appears  the  mountain 
range  may  have  a  stabilizing  effect  on  the  growing  distur¬ 
bance.  The  required  25  mb  pressure  difference  is  reached  at 
84  h.  Since  a  twelve  hour  period  is  required  to  create  the 
valley,  growth  will  be  initiated  at  72  a  model  time  in  sub¬ 
sequent  runs.  Since  the  valley  must  be  located  east  of  the 
anticyclone  in  order  to  taka  advantage  of  the  associated 
northwest,  flow,  the  western  edge  of  tae  valley  is  placed 
directly  south  of  the  leading  portion  of  northwesterly  winds 
at  84  h. 

Development  of  the  disturbance  in  the  potential  tempera¬ 
ture  and  wind  fields  is,.,  very  similar  to  the  first  run. 
There  is  a  maximum  in  potential  temperature  over  the  moun¬ 
tains  but  this  is  simply  due  to  the  increased  elevation. 

The  small  scale  noise  apparent  in  the  potential  tempera¬ 
ture  fields  of  the  first  run  also  appear  in  this  experiment 
although  at  a  later  time,  possibly  due  to  the  much  slower 
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development  of  the  cyclone-anticyclone  system.  The  strong 
cyclogenesis  and  intense  pressure  gradients  also  appear 
later  in  the  second  model  run.  Wave  development  is  approxi¬ 
mately  six  hours  slower  when  the  mountains  are  included  and 
this  corresponds  to  the  time  lag  observed  in  the  cyclcgene- 
sis  and  noise  growth. 

It  should  be  noted  that  potential  temperature  is  used  in 
the  analysis  of  model  results  since  its  fields  will  not 
include  effects  of  the  adiabatic  processes  occurring  in  flow 
over  mountains. 

B.  EXPERIMENT  3 

The  initialization  for  this  experiment,  and  all  subse¬ 
quent  experiments,  is  provided  by  the  same  analytic  scheme 
used  in  the  control  runs.  Terrain  consists  of  an  east-west 
mountain  range  centered  at  30°N.  The  valley,  built  between 
72  and  84  h,  is  twelve  gridlangths  wide  or  nearly  one-third 
the  length  of  the  mountain  range  (Fig.  12)  .  The  critical 
time  period  for  the  occurrence  of  a  surge  will  be  the  first 
24  h  after  the  valley  is  completed  (84-133  h) .  It  is  during 
this  period  that  the  anticyclone  is  located  northwest  of  the 
valley. 
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Development  of  the  disturbance  through  72  h  is  described 


by  control  run  *2.  The  sea-level  pressure  field  at  84  h 
shows  a  trough  extending  across  the  valley  and  along  the 
northern  slope  of  the  eastern  range.  A t  the  same  time  a 
small,  high  pressure  ridge  can  be  seen  in  the  northwest 
region  of  the  valley.  For  the  next  twelve  hours  the  ridge 
continues  to  build  southward  until,  at  96  a,  it  extends  the 
length  of  the  western  valley  wall  splitting  the  trough  into 
two  separate  ]ow  centers.  The  pressure  gradient  in  the 
tropical  region  south  of  the  valley  has  also  increased  (Fig. 
13) .  As  the  ridge  continues  to  expand  southward,  the  east¬ 
ern  portion  of  the  trough  develops  a  very  strong,  tight 
pressure  gradient  in  the  valley  and  along  tae  northern  moun¬ 
tain  slopes  (Fig.  14).  The  ridge  continues  to  strengthen 
and  remains  located  along  the  western  wail  through  114  h. 
After  this  time  the  anticyclone  moves  eastward  to  a  position 
due  north  of  the  valley  and  the  ridge  begins  to  weaken. 

Also  noted  in  the  results  is  an  oscillation  in  the 
intensity  of  the  pressure  field.  It  begins  when  valley 
growth  is  initiated  and  has  a  period  of  approximately  24  h. 
The  observed  oscillation  is  possibly  the  result  of  gravity 
waves  set  up  by  the  relatively  sudden  (12  h)  creation  of  the 
valley. 
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The  second  stage  of  the  monsoon  surge  was  earlier 
described  as  being  characterized  by  a  sharp  drop  in  surface 
dewpoint.  Since  this  is  a  dry  model,  and  as  observational 
studies  (Chang  et  al,  1979,  Chu,  1978)  show  this  stage  can 
also  be  associated  with  a  drop  in  surface  temperatures,  the 
potential  temperature  fields  are  analyzed  for  possible  surge 
characteristics.  (Beyond  this  point  'temperature'  refers  to 
potential  temperature  unless  otherwise  stated.) 

In  agreement  with  the  development  of  the  surface  pres¬ 
sure  ridge,  a  weak  tongue  of  cold  air  has  formed  along  the 
western  valley  wall  at  84  h.  However,  no  similar  gradient 
exists  along  the  astern  wall  (Fig.  15). 

Twelve  hours  after  the  valley  is  completed,  tne  cold  air 
has  reached  the  southern  valley  boundary.  Minimum  tempera¬ 
ture  at  24°H  has  dropped  from  310K  at  34  h  to  303K  at  96  h. 
A  tight  temperature  gradient  exists  along  the  western  side 
of  the  cold  air  intrusion  and  also  along  the  northern  moun¬ 
tain  slopes  where  temperatures  have  remained  relatively  con¬ 
stant.  The  cold  air  tongue  can  easily  by  seen  at  the  two 
lowest  levels  (Fig.  16)  .  At  the  mid-levels  a  very  slight 
tightening  of  the  gradient  can  be  seen  over  the  western  val¬ 
ley  wail  when  compared  to  the  eastern  wall,  but  only  over 
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the  northern  half  of  the  valley.  At  the  highest  two  levels, 
gradients  over  the  two  valley  walls  are  nearly  identical  so 
that  the  apparent  changes  in  temperature  are  probably  due  to 
the  effect  of  terrain  variation  on  the  slope  cf  the  sigma 
surface. 

The  cold  air  continues  to  penetrate  southward  for  the 
next  twelve  hours  reaching  its  maximum  southward  extent  at 
138  h  (Fig.  17).  At  24° N  the  minimum  temperature  is  298K,  a 
-12K  change  from  84  h.  The  cold  air  does  not  penetrate  much 
beyond  the  southern  valley  border.  At  23°N  the  temperature 
change  is  approximately  -2K  in  24  h.  While  southward  move¬ 
ment  of  the  cold  air  ceases  at  108  h,  it  does  expand  east¬ 
ward,  apparently  in  response  to  the  eastward  movement  of  the 
anticyclone  located  north  of  the  valley.  The  cold  air 
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At  the  higher  levels,  sigma  =  7/12  and  3/12,  a  small 
tongue  of  cold  air  forms  between  78  and  34  a  above  the  west¬ 
ern  valley  wall.  The  southward  extent  is  slightly  less  than 
that  at  the  lowest  level.  As  the  low-level  cold  air  tongue 
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moves  rapidly  southward  after  84  h,  the  cold  air  at  these 
levels  moves  very  slowly,  only  four  degrees  between  34  and 
102  h  at  sigma  =  7/12.  Little,  if  any,  eastward  expansion 
is  seen.  Beyond  108  h  the  very  slight  gradients  are  diffi¬ 
cult  to  discern,  primarily  due  to  the  increasing  small  scale 
noise.  At  sigma  =  9/12  the  cold  air  more  closely  follows 
the  behavior  of  the  cold  air  tongue  at  the  lowest  level 
showing  continued  south  and  eastward  expansion  on  the  same 
time  scale,  although  of  less  intensity. 

To  verify  that  the  changes  in  the  potential  temperature 
fields  at  the  upper-levels  are  essentially  due  to  the  use  of 
the  sigma  coordinate  system,  the  temperature  (I)  fields  were 
analyzed  at  350,  700  and  500  mb  between  8u  and  102  h.  At 
850  mb  the  intrusion  cf  cold  air  can  easily  be  traced.  At 
700  mb  the  cold  air  tongue  is  not  discecnable  although  this 
may  partially  be  due  to  the  adiabatic  warming  effects  appar¬ 
ent  at  both  the  700  and  850  mb  levels.  At  500  mb  some  small 
scale  noise  is  introduced  by  flow  over  and  around  the  valley 
but  there  is  no  evidence  of  the  southward  movement  cf  cold 
air.  Therefore  it  is  reasonable  to  assume  that  the  changes 
in  the  potential  temperature  fields  at  the  higher  levels  are 
due  to  the  effect  terrain  variation  has  on  the  slope  of  the 
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flows.  The  strong,  cross-isobaric  southwest  winds  in  the 
eastern  valley  coincide  with  a  very  tignt  pressure  gradient 
and  may  represent  a  sloping  front  in  the  valley. 

At  96  h  the  north  wind  regime  extends  the  length  of  the 
western  valley  wall  and  is  spreading  eastward  across  the 
valley  (Fig.  19).  At  102  h  the  north  winds  have  reached 
their  maximum  extent;  the  wind  regime  occupies  fully  half 
the  valley  but  dees  net  extend  south  of  it.  The  flow  along 
the  valley  wall  is  no  longer  age ostrophic.  The  north  winds 
entering  the  valley  between  84  and  102  a  are  identifiable 
only  at  the  very  lowest  levels,  in  agreement  with  the  defi¬ 
nition  of  a  monsoon  surge  as  a  low-level  phenomenon.  By  103 
h  the  anticyclone  has  moved  sufficiently  eastward  that  the 
light  northerly  winds  associated  with  tee  ridge  line  have 
replaced  the  strong  northwest  flow  previously  located  north 
of  the  valley  (Fig.  20).  The  flew  into  the  valley  is 
greatly  reduced  and,  without  this  support,  the  surge  weakens 
and  dissipates.  The  eastern  region  continues  to  be  domi¬ 
nated  by  southwest  flow. 

It  is  clear  that  the  movement  of  the  north  winds  through 
the  valley  correlates  well  with  the  development  of  the  cold 
air  tongue  and  pressure  ridge  discussed  earlier. 
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In  summary,  saveral  observations  can  be  made  from  the 
analysis  of  these  fields.  First,  all  three  fields  shew  the 
vary  rapid  southward  movement  of  a  cold,  high  pressure  air- 
mass  along  the  western  valley  wall  batwear.  84  and  96  h. 
Much  slower  southward  movement  continues  through  108  h. 
Second,  the  airmass  expands  eastward  across  the  valley  in 
response  to  the  continued  eastward  movement  of  the  anticy¬ 
clone.  And  third,  the  airmass  never  penetrates  the  tropical 
regions  but  remains  within  the  confines  of  the  valley. 

C.  EXPERIMENT  4 

In  this  experiment  the  mountain  range  is  again  centered 
at  3  0°N .  The  valley,  built  between  72  and  84  h,  is  24 
grid-lengths  wide,  ie.  two-thirds  of  me  mountain  range  in 
length.  The  terrain  then  is  actually  the  small  east-west 
oriented  mountain  range  displayed  in  Fig.  21.  As  in  Experi¬ 
ment  #3,  the  critical  period  for  the  surge  is  between  82  and 
108  h  while  the  anticyclone  is  located  northwest  of  the 
valley. 

The  first  72  h  of  development  are  described  by  Control 
Run  #2.  As  in  the  previous  experiment  the  creation  of  the 
valley  causes  an  oscillation  in  the  pressure  (ie.  mass) 
field.  The  oscillation  in  this  experiment  is  larger  than 
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-hat  of  the  previous  case  since  the  reduction  of  certain  is 


greater.  The  rapid  increase  in  cyclogenesis  also  occurs 
earlier,  supporting  the  assumption  that  the  mountains  exert 
a  stabilizing  influence  on  the  disturbance. 

At  84  h  a  small  southward  extension  of  the  hign  pressure 
ridge  can  be  seen  over  the  eastern  mountain  slope  (western 
valley  wall).  A  low  is  centered  over  the  northern  valley 
with  a  weak  trough  extending  southwest  over  the  valley  wall. 
A  tight  gradient  develops  south  of  the  low,  extending  over 
the  eastern  valley  wall  and  continuing  along  the  northern 
mountain  slopes.  The  trough  is  strongest  at  78  h  and  is 
subsequently  pushed  southward  by  the  expanding  high  pressure 
ridge.  The  ridge  advances  southward  during  the  next  twelve 
hours  reaching  the  southern  valley  edge  at  96  h  (Fig.  22)  . 
Thus  far  movement  of  the  ridge  is  quite  similar  to  Experi¬ 
ment  #3,  although  it  does  appear  to  be  slightly  weaker  and 
not  as  readily  discernaole. 

Between  96  and  108  h  the  mid-latitude  ridge  weakens  con¬ 
siderably.  Concurrently,  a  strong  pressure  gradient  builds 
along  the  western  edge  of  the  ridge,  in  part  due  to  a  trough 
extending  southward  from  the  intense  cyclone  located  to  the 
northwest  (Fig.  23).  A  similar  gradient  is  seen  during  the 
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same  period  in  Experiment  #3  but  is  nor  nearly  as  intense. 
At  both  102  and  108  h  there  is  evidence  that  the  ridge 
extends  slightly  beyond  the  edge  of  the  valley. 

The  development  and  movement  of  the  ridge  agrees  very 
closely  with  that  of  the  earlier  experiment  using  a  smaller 
valley.  In  both  cases  the  ridge  begins  development  while 
the  valley  is  growing,  moves  southward  along  the  western 
wail  at  the  same  rate,  reaches  the  southern  edge  at  the  96  h 
point  and  weakens  when  the  anticyclone  moves  to  a  location 
due  north  of  the  valley.  There  is,  however,  a  slight  dif¬ 
ference  in  intensity.  The  third  experiment  shows  a  well 
defined,  easily  discernable  ridge  from  34  to  108  h.  In  the 
present  experiment  the  gradient  along  the  western  side  of 
the  ridge  is  greater,  but  the  ridge  itself  is  not  as  well 
defined. 

A  weak  tongue  of  cold  air  can  be  seen  along  the  western 
valley  wall  at  84  h  (Fig.  24)  .  An  increase  an  the  tempera¬ 
ture  gradient  is  apparent  in  the  extreme  western  valley  at 
the  lower  three  levels  although  it  is  very  weak  at  the 
uppermost  of  these  (sigma  =  7/12).  These  levels  correspond 
to  pressure  heights  lower  than  600  mbs  in  a  standard  atmcs- 
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Through  the  next  twexve  hours  the  cold  air  tongue 


moves  southward  and  strengthens.  By  96  h  the  cold  air  has 
reached  the  southern  edge  of  the  valley  (Fig.  25).  The  cold 
air  tongue  is  much  better  defined  than  the  pressure  ridge. 
It  is  nearly  identical  to  the  cold  air  tongue  seen  in  the 
western  valley  for  Experiment  #3.  The  minimum  reaperature 
at  24°N  has  fallen  from  3  10K  at  84  h  to  300K  at  96  h.  The 
temperature  change  is  primarily  limited  to  the  western  third 

of  the  valley,  temperatures  in  the  eastern  two-thirds  have 

« 

remained  relatively  constant  following  completion  of  the 
valley. 

The  cold  air  continues  to  move  south  and  eastward 
through  the  next  twelve  hours.  Maximum  southward  extant  is 
reached  at  138  h.  At  24°n  the  minimum  temperature  is  298K 
and  represents  the  same  change  in  temperature,  -12K,  as  was 
seen  in  Experiment  *3  (Fig.  26)  .  While  southward  movement 
of  the  cold  air  has  ceased,  if  does  continue  to  expand  east¬ 
ward  well  past  this  time.  The  primary  axis  of  the  cold  air 
tongue  remains  along  the  western  valley  wall  and  retains  its 
identity  through  120  h.  A  pool  of  cool  air  remains  in  the 
western  valley  for  another  twelve  hours.  At  sigma  =  9/12 
the  cold  air  is  discernable  wail  past  108  h.  At  nigher  lev¬ 
els  the  small  scale  noise  has  more  effect  and  the  cold  air 
is  difficult  to  trace  past  this  time. 
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The  upper-levels  display  the  same  behavior  seen  in 
Experiment  #3.  A  slight  increase  in  eastward  expansion  at 
these  levels  is  seen,  possibly  due  to  the  much  wider  valley. 

As  with  the  pressure  and  temperature  fields  the  develop¬ 
ment  of  the  north  wind  surge  in  this  experiment  closely  fol¬ 
lows  that  of  Experiment  #3.  At  84  h  the  north  winds  have 
reached  mid-valley  and  continue  to  spread  rapidly  southward. 
Again  the  leading  edge  of  the  wind  surge  is  a geost rcphic.  A 
strong  southwesterly  flow  has  set  up  in  the  valley  appar¬ 
ently  due  to  the  tight  pressure  gradient  located  there  (Fig. 
27) .  The  small  cyclones  seen  in  Experiment  #3  are  not 
apparent  here  although  cyclonic  flew  exists,  possibly  due  to 
the  convergence  of  the  southwest  flow  from  the  tropics  and 
the  opposing  flow  out  of  the  mid-latitude  anticylones.  This 
cyclonic  flow  weakens  considerably  after  34  h. 

At  96  h  the  north  wind  regime  extends  almost  the  entire 
length  of  the  valley  and  is  spreading  eastward.  The  south¬ 
ward  extent  of  the  north  winds  along  the  valley  wall  is 
identical  to  that  of  Experiemtn  *3  (Fig.  28).  Further  east 
the  perimeter  of  the  north  wind  regime  is  nearly  two  degrees 
farther  south  of  that  reached  in  the  earlier  experiment. 
The  mid-latitude  ridge  weakens  during  the  next  srx  hours. 
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and  with  its  continued  eastward  movement, 
decrease  in  the  northwesterly  winds  supporting  the  surge. 
Thus  only  a  very  little  southward  movement  takes  place 
between  96  and  102  h,  the  time  of  maximum  southward  extent 
of  the  north  winds.  The  ageostrophic  component  of  the  flow 
weakens  rapidly  as  southward  movement  along  the  valley  wall 
slows.  The  winds  become  much  weaker  but  do  not  decrease  in 
coverage  through  108  h  (Fig.  29).  A  slight  eastward  move¬ 
ment  is  seen  as  the  ridge  moves  past  the  mountains  and  north 
of  the  valley.  After  108  h  the  north  winds  diminish  and 
move  eastward  in  response  to  the  weakening  and  continued 
eastward  movement  of  the  ridge.  After  120  h  the  wind  surge 
is  no  longer  discernable. 

The  surge  is  identifiable  only  at  tne  two  lowest  levels. 
At  sigma  =  9/12  the  surge  is  weaker  and  less  extensive  than 
at  the  lowest  level.  It  dissipates  very  rapidly  after  102  h 
as  the  ridge  at  this  level  moves  north  of  the  valley. 

In  this  experiment  a  slight  increase  in  the  eastward 
extent  of  the  north  winds  is  noted  when  compared  to  Experi¬ 
ment  #3.  Since  topography  is  the  only  parameter  changed,  it 
is  reasonable  to  assume  the  increase  in  eastward  extent  is  a 
result  of  the  expanded  valley.  It  should  be  noted  though 
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that  while  the  valley  size  has  doubled,  the  corresponding 
increase  in  longitudinal  coverage  by  the  north  winds  is 
extremely  small. 

Several  observations  can  be  made  from  the  results  of 
Experiment  #4.  First,  all  three  fields  show  the  rapid 
southward  movement  of  the  cold,  high  pressure  airmass  along 
the  western  valley  wall  between  34  and  96  h.  Subsequent 
southward  expansion  is  much  slower.  Second,  the  airmass 
expands  eastward  across  the  valley  in  response  to  the  east¬ 
ward  movement  of  the  anticyclone.  And  third,  there  is  a 
tendancy  towards  cooling  and  rising  pressures  just  south  of 
the  valley  although  the  ncrth  winds  remain  within  the  valley 
confines. 


D.  EXPERIMENT  5 

In  this  experiment  the  mountain  range  is  centered  at 
33  °N ,  eight  degrees  farther  north  than  in  previous  experi¬ 
ments.  The  valley  is  twelve  gridlsngths  wide,  as  in  Experi¬ 
ment  #3.  Terrain  is  shewn  in  Fig.  33.  Cyclogenesis  is 
considerably  slower  than  that  seer,  in  either  of  Experiments 
3  or  4.  This  is  due  to  the  placement  of  the  northern  moun¬ 
tain  slopes  at  the  same  latitude  as  the  maximum  amplitude  of 
the  disturbance.  Srowth  cf  the  perturbation  is  prevented  by 
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affectively  blocking  the  north-south  heat  flux.  For  this 
reason  the  initiation  of  valley  growth  was  ielayed  six  hours 
to  allow  additional  time  for  development  of  the  disturbance. 
The  growth  period  for  the  valley  is  now  between  78  and  90  h. 
As  with  the  previous  experiments  the  creation  of  the  valley 
causes  an  oscillation  in  the  pressure  field,  the  approximate 
amplitude  of  that  in  Experiment  #3. 

At  90  h  a  strong  low  is  located  in  the  northeast  of  the 
valley  trough  extending  southwest  over  the  valley  wall. 
This  results  in  a  tight  pressure  gradient  in  the  eastern 
valley  and  along  the  northern  mountain  slopes  similar  to  the 
gradients  found  in  previous  experiments.  The  mid-latitude 
anticyclone  is  neither  as  strong  nor  as  well  developed  as  in 
Experiments  3  and  4,  and  although  there  does  appear  to  be  a 
slight  southward  extension  of  the  high  pressure  along  the 
western  valley  wall,  it  is  very  weak.  During  the  next 
twelve  hours  a  ridge  of  high  pressure  builds  along  the  val¬ 
ley  wail,  however  it  is  difficult  to  determine  whether  the 
high  pressure  originates  in  the  anticyclone  to  the  northwest 
or  in  the  tropical  regions  south  of  the  valley  (Fig.  31). 
As  the  analysis  of  the  wind  fields  will  later  show,  the  high 
pressure  air  does  originate  in  the  mid-latitude  anticyclone. 
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At  102  h  the  anticyclone  is  ill-defined  and  the  mountain 
region  to  its  south  populated  by  small  high  and  low  centers. 

Between  102  and  114  h  the  ridge  moves  east  into  the  val¬ 
ley  and  strengthens,  this  is  particularily  noticeable  at  114 
h  (Pig.  32).  It  appears  the  anticyclone  has  slipped 
beneath  the  cyclone  and  is  re-forming  in  the  northern  val¬ 
ley.  During  this  time  a  weak  trough  can  be  seen  forming 
over  the  northwestern  valley  wall. 

At  120  h  the  anticyclone  is  located  in  the  northeast 
valley  region  with  a  weak  ridge  extending  southwest  througn 
the  valley  and  a  stronger,  better  defined  ridge  to  the 
northwest.  Low  pressures  are  found  over  the  western  valley 
wall  due  to  a  weak  trough  forming  ahead  of  the  mid-latitude 
cyclone. 

In  earlier  experiments  the  high  pressure  ridge  developed 
and  strengthened  along  the  western  valley  wall  only  to  dis¬ 
sipate  as  the  anticyclone  moved  to  a  position  northeast  of 
the  ridge.  In  this  experiment  the  ridge  is  weaker,  net 
well-defined  and  moves  eastward  with  the  anticyclone.  Since 
the  northern  valley  is  located  at  the  latitude  of  maximum 
amplitude  for  the  disturbance,  creating  the  valley  allows 
the  anticyclone  to  reform  in  this  area  as  the  weak  ridge 
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moves  eastward  into  the  valley.  The  same  reasoning  can  be 
used  to  explain  the  strong  development  experienced  by  the 
low  during  the  valley  growth  period;  it  was  located  along 
the  northern  mountain  slopes  at  the  same  longitude  as  the 
valley. 

A  weak  tongue  of  cold  air  can  be  seen  along  the  western 

valley  wall  at  90  h  (Fig.  33).  This  gradient  is  found  in 

the  lowest  three  levels,  weakening  with  height.  At  the 
upper-levels  (sigma  =  5/12,  3/12,  1/12)  nearly  identical 

gradients  exist  at  both  ends  of  the  valley  and  are  due  to 
the  effects  of  terrain  on  the  slope  of  sigma  surfaces.  The 
cold  air  tongue  strengthens  and  moves  rapidly  southward  in 
the  following  six  hours.  3y  10  2  h  the  cold  air  has  reached 

the  southern  edge  of  the  valley  (Fig.  34)  .  The  development 

and  movement  of  the  cold  air  tongue  is  much  better  defined 
than  that  cf  the  pressure  ridge.  It  is  nearly  identical  in 
development  and  intensity  to  the  cold  air  tongues  seen  in 
Experiments  3  and  4.  The  minimum  temperature  at  32°N  (a 
comparable  latitude  to  24°N  relative  to  the  position  of  the 
mountains)  is  295K  compared  to  306K  at  90  h.  This  is  a 
temperature  change  of  —  1  1 K  compared  to  -12K  and  -10K  for 
Experiments  3  and  4  respectively.  Again  the  temperature 


change  is  primarily  limited  to  the  western  third  of  the 
valley. 

The  cold  air  continues  to  move  south  and  eastward  during 
the  next  twelve  hours,  reaching  its  maximum  southward  extent 
at  1  14  h  (Fig.  35)  -  At  32°N  the  minimum  temperature  is 
293K,  a  temperature  change  of  —  1 3 K  compared  to  -12K  in 
Experiments  3  and  4  for  a  24-h  period  following  completion 
of  the  valley.  Eastward  expansion  of  the  cold  air  continues 
past  114  h.  In  the  previous  experiments  the  axis  of  the 
cold  air  tongue  remained  along  the  western  valley  wall.  In 
this  experiment  the  cold  air  moves  eastward  as  well,  pacti¬ 
on  la ril y  in  the  northwest  region.  This  is  probably  the 
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similar  movement  is  seen  at  sigma  *  7/12.  The  increased 
small  scale  noise  at  these  levels  mattes  analysis  of  the 
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temperature  fields  difficult  past  114  h. 


The  highest  levels 


(sigma  =  5/12  to  1/12)  show  no  sigr.  of  southward  movement  of 
cold  air  at  any  time. 

As  in  the  previous  experiments  southward  deflection  of 
the  northwest  winds  into  the  valley  occurs  before  the  valley 
is  completed.  3y  90  h  the  north  winds  have  reached  aid-val¬ 
ley  with  a  shear  line  developing  over  the  valley  floor  in 
advance  of  the  northerly  flow.  North  winds  along  the  valley 
wall  are  strongly  ageostrophic.  Again  a  strong  southwest 
flow  has  set  up  in  response  to  the  tigat  pressure  gradient 
in  the  eastern  valley.  Although  the  valley  is  the  same  as 
that  used  in  Experiment  #3,  the  subsequent  pattern  of  devel¬ 
opment  is  slightly  different.  At  96  h  the  northwest  winds 
reach  'the  eastern  valley  wall,  whereas,  in  Experiment  *3  the 
strong  southwest  flow  in  this  region  prevents  the  northwest 
winds  from  spanning  the  valley.  The  wind  surge  is  actually 


strongest  at  this 

ti 

me,  and  has 

its  greatest 

southward 

extent  (Fig.  36)  . 

By 

102  h  the  ridge 

is  north  of 

the  valley 

wall  and  the  winds 

are 

beginning  to 

weaken,  the 

flow  is  no 

longer  ageostrophic 

• 

The  north  winds 

never  reach 

the  south- 

ern  edge  of  the  valley 

(Fig.  37). 

The  surge  continues  to 

weaken  rapidly  and 

at 

108  h  the  anticyclone  is 

reforming 
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within  the  valley.  North  winds  axis-1:  in  tne  central  valley 
region  but  dc  not  coincide  with  the  positon  of  tne  cold, 
high  pressure  ridge.  Six  hours  later  the  anticyclone  is 
established  over  the  north  central  valley  and  no  evidence  of 
the  surge  exists  (Fig.  3  8).  It  is  clear  that  positioning 
tne  mountain  range  at  38°N  has  a  detrimental  effect  on  the 
development  of  the  surge. 

The  surge  is  identifiable  at  only  the  two  lowest  levels. 
At  sigma  ~  9/12  the  surge  is  less  extensive  and  weaker  than 
at  the  lowest  level.  At  sigma  =  7/12  and  5/12  a  weak 
cyclonic  turning  of  the  winds  is  apparent  but  r.o  definite 
southward  flow  can  be  seen.  The  cyclonic  flow  disappears  as 
the  anticyclone  reforms  over  tne  valley. 

In  summary  the  following  observations  can  be  made  from 
this  experiment.  First,  all  three  fields  show  the  rapid 
southward  movement  of  cold,  high  pressure  air  along  the 
western  valley  wall  during  the  first  six  hours  following 
valley  completion.  Second,  the  movement  and  development  of 
the  cold  air  tongue  is  much  better  defined  than  that  of  the 
pressure  ridge  and  wind  surge,  both  of  which  weaken  consid¬ 
erably  after  102  h.  Third,  the  pressure  ridge  and  wind 

the 


surge  axis  move  eastward  off 


he  valley  wail  due  to 


development  of  the  anticyclone  within  tne  valley,  only  the 
northern  most  portion  of  the  cold  air  tongue  moves  eastward. 
And  fourth,  although  there  is  a  slight  oool;ng  and  pressure 
rise  south  of  the  valley,  the  north  winds  remain  well  within 
the  valley. 

S.  EXPERIMENTS  6  AND  7 

The  valley  was  included  as  part  of  the  initial  terrain 
configuration  for  these  two  experiments.  The  mountains  are 
centered  at  30°N  and  the  valley  is  eignt  gridlengths  wide, 
approximately  one-quarter  the  length  of  the  mountain  range. 
This  is  the  original  mountain -valley  configuration  used  with 
the  model,  the  larger  valley  was  used  m  subsequent  runs  to 
avoid  interference  between  the  strong  southwesterly  flow  in 
the  eastern  valley  and  developing  soutaward  flow  in  the 
west.  The  southward  movement  of  a  cold,  high  pressure  air- 
mass  occurs  with  both  valley  sizes,  however  the  airmass  is 
more  easily  seen,  and  so  more  easily  analyzed,  with  the 
larger  valley. 
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As  expected  no  "surge"  occurs.  Tde  current  that  dees 
form  in  the  valley  is  southwesterly,  flowing  from  high  to 
low  pressures.  The  strong  southwesterly  flow  seen  in  Exper¬ 
iments  3  through  5  is  in  part  due  to  the  mean  pressure  gra¬ 
dient  that  exists  between  the  northern  and  southern  valley 
borders.  3y  30  h  a  small  cyclone  has  formed  over  the  west¬ 
ern  valley  wall,  probably  due  to  the  conservation  of  poten¬ 
tial  vorticity.  This  indicates  that  the  small  cyclone  seen 
here  and  in  ether  experiments  is  a  result  of  the  mean  flow 
interacting  with  the  terrain. 

Experiment  #7  is  initialized  with  the  perturbation. 
During  the  first  portion  of  the  experiment,  while  the  dis¬ 
turbance  is  weak,  southwesterly  flow  forms  in  the  valley  in 
response  to  the  existing  north-south  pressure  gradient. 
Again  at  30  h  a  small  cyclone  has  formed  over  the  western 
valley  wall.  During  the  next  twelve  hours  the  center  of  the 
cyclone  moves  ever  the  valley  floor  althougn  the  north  winds 
remain  located  over  the  valley  wall.  The  mid-latitude  dis¬ 
turbance  continues  to  strengthen  and  by  66  h  the  northerly 
winds  are  being  enhanced  by  the  northwest  flow  ahead  of  the 
approaching  ridge.  The  northerly  flow  strengthens  as  the 
ridge  moves  eastward,  at  96  h  it  extends  slightly  beyond  the 
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valley  edge,  however  this  is  its  maximum  southward  extent. 
Subsequently  the  mid-latitude  ridge  axis  moves  north  of  the 
valley  and  the  supporting  northwest  flow  diminishes.  The 
flow  along  the  valley  wall  is  never  ageostrophic  during  this 
time. 

The  cyclone  is  readily  apparent  in  the  sea-level  pres¬ 
sure  fields.  However,  not  until  the  ridge  approaches  the 
valley  does  a  high  pressure  ridge  build  along  the  valley 
wall  (72  h)  .  The  ridge  weakens  as  the  anticyclone  moves  tc 
a  position  north  of  the  valley. 

The  potential  temperature  field  shows  very  little  cool¬ 
ing  along  the  western  valley  wall  compared  to  previous 
experiments.  minimum  temperature  along  24°H  at  90  h  is  299K 
compared  to  303K  at  66  h,  a  decrease  of  only  in  2a  h.  A 
very  small  temperature  change  when  compared  to  Experiments  3 
through  5. 
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gure  22.  Sea-Level  Pressure,  3u  and  96  h 
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The  interaction  of  a  aid-latitude  baroclinic  wave  with 
several  different  mountain- valley  configurations  was  exam¬ 
ined.  The  purpose  was  to  determine  the  contribution  made  by 
a  well-developed  baroclinic  wave  to  the  initiation  of  mon¬ 
soon  surges.  It  was  found  that  these  waves  can  initiate  a 
monscon  surge  although  such  surges  are  wea<c  and  limited  in 
extent. 

The  mountains  in  each  experiment  are  an  east-west  range 
located  south  of  the  mid-latitude  wave.  As  was  first  noted 
in  Experiment  #2  the  mountains  appear  to  have  a  stabilizing 
affect  on  the  growth  of  the  wave.  The  effect  is  especially 
evident  when  the  mountains  are  placed  farther  north,  nearer 
the  center  of  the  disturbance.  The  reduced  growth  rate  is 
probably  caused  by  the  blocking  of  the  v-component  of  the 
wind  by  the  mountain  barrier  to  the  sou‘h.  This  would  slow 
the  wave  growth  because  the  energy  conversion  is  porpor- 
tional  to  V'9'.  Walker  (  1982)  also  found  this  effect  with  a 
disturbance  moving  along  the  north  side  of  a  mountain  range. 

Previously,  several  observations  were  made  about  the 
nature  of  monscon  surges.  These  include:  (1)  the  occurrence 
of  rapid,  ageostrophic  flow  from  high  to  low  pressure,  (2) 
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strengthening  of  the  north-northeasterly  winds,  (3)  a  sharp 
drop  in  surface  temperature  and  (4)  the  shallow  structure  of 
the  surge.  Tc  seme  extent,  each  of  these  characteristics  is 
simulated  by  cne  or  more  of  the  experiments. 

Experiment  #6,  in  which  only  the  zonal  mean  flow  inter¬ 
acts  with  the  terrain,  shews  no  evidence  of  a  surge.  All 
other  experiments,  in  which  the  perturbation  is  added  to  the 
mean  flow,  do  produce  weak  monsoon  surges.  This  indicates 
the  eastward- mov ing  anticyclone  associated  with  the  mid-la¬ 
titude  baroclinic  wave  is  the  source  of  the  surges  occurring 
in  this  study. 

In  Experiments  3  and  4  the  mountains  are  centered  at 
30  °N ,  the  western  edge  of  each  valley  is  located  at  the  sane 
point  although  the  fourth  valley  is  twice  the  size  of  the 
third.  In  each  case  a  cold,  high  pressure  airmass  moves 
southward  along  the  western  valley  wall  during  the  twelve 
hours  following  completion  of  the  valley.  The  leading  edge 
of  the  surge  is  ageostrophic. 

Figs.  39  and  40  show  the  change  of  potential  temperature 
with  time  in  the  valley  for  Experiments  3  and  4.  The  pro¬ 
files  extend  from  14°n  to  46° N  along  the  western  valley 
floor  (55°  long),  the  valley  is  located  between  22°N  and 
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38°N.  In  each  case  little,  if  any,  change  in  temperature 
occurs  at  14°N  although  at  18°N  the  temperature  decreases  7 K 
between  96  and  108  h.  In  the  valley  the  maximum  temperature 
change  occurs  earlier,  between  84  and  96  h,  and  is  much 
larger,  -16K  and  -17K  at  26°N.  The  surge  is  well  defined 
within  the  valley  but  dissipates  rapidly  south  of  22°N. 
This  is  indicative  of  the  we  a  ic  character  of  the  monsoon 
surges  in  this  study  since  the  cold  air  in  observed  cases 
has  been  traced  well  into  the  tropical  regions  while  in 
these  experiments  it  cannot  be  found  even  eight  degrees 
south  of  the  valley. 

As  can  be  seen  the  results  of  the  two  experiments  are 
nearly  identical,  indicating  the  size  of  the  valley  is  not 
an  important  consideration. 

-  oscillation  in  the  pressure  (mass)  field  is  appar¬ 
ently  a  gravity  wave,  the  result  of  an  inward  rush  of  air 
filling  the  "vacuum"  produced  by  the  creation  of  the  valley. 
It  has  roughly  a  24-h  period,  which  may  be  related  to  the 
twelve  hour  period  in  which  the  valley  is  built,  and  it 
slowly  dampens  during  the  subsequent  60  h  of  the  experiment. 
It  has  no  observable  effect  on  the  pressure  gradients  and 
was  not  evident  in  either  the  potential  temperature  or  the 
wind  fields.  The  oscillation  also  occurs  in  Experiment  #5. 
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In  Experiment  #5  the  rapid  southward  movement  of  the 
cold  airmass  occurs  only  during  the  first  six  hours  after 
valley  completion.  This  is  due  to  the  placement  of  the 
mountains  at  38°N  resulting  in  a  much  weaker  disturbance. 
Fig.  41  shows  the  change  in  potential  temperature  with  time 
in  the  valley.  The  valley  is  located  between  30°N  and  46°N, 
the  profile  is  along  the  western  valley  floor. 

The  decrease  in  temperature  is  evident  in  the  valley 
between  90  and  102  h.  As  in  Experiments  3  and  4  the  cold 
air  can  be  traced  only  four  degrees  south  of  the  valley. 
Unlike  the  two  previous  experiments  warming  occurs  during 
the  second  twelve  hour  period,  due  to  the  re-development  of 
the  anticyclone  within  the  valley.  The  southerly  winds  in 
the  western  region  of  the  anticyclone  prevent  the  flow  of 
cold  air  into  the  valley  and  effectively  "shut  off"  the 
surge. 

The  mcnsoor.  surges  initiated  in  these  experiments  are 
evident  only  in  the  lowest  sigma  levels,  at  sigma  ~  11/12, 
9/12  and  7/12.  In  a  standard  atmosphere,  over  flat  terrain, 
these  correspond  to  pressure  heights  less  than  675  mb.  The 
surge  is  strongest  at  the  lowest  level  and  dissipates  with 
height,  becoming  very  weak  at  sigma  =  7/12.  This  is 
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observed  in  both  the  wind  and  potential  temperature  fields 
in  all  experiments  and  agrees  closely  with  the  definition  of 
the  monsoon  surge  as  a  low-level,  or  surface,  phenomenon. 

These  results  suggest  that  the  aid-latitude  baroclinic 
wave  can  initiate  monsoon  surges.  The  surges  are  weak  and 
limited  to  the  valley  but  otherwise  exhibit  the  characteris¬ 
tics  of  observed  monsoon  surges.  The  weakness  of  these 
surges  indicates  other  forcing  mechanisms  are  necessary  in 
order  to  adequately  model  the  monsoon  surge.  These  may 
include  a  strong,  stationary  high  located  northwest  of  the 
valley  similar  to  the  anticyclone  associated  with  the  winter 
monscon.  other  analytical  studies  indicate  that  an  upper- 
tropospheric,  planetary  scale  wave  number  three  is  necessary 
for  the  occurrence  of  strcng  monsoon  surges.  Such  a  wave 
positioned  with  the  ridge  to  the  west  of  the  surface  anticy¬ 
clone  would  result  in  strong  northwesterly  flow  over  the 
high  pressure  center.  This  would  initiate  cold  air  advec- 
tion  throughout  the  troposphere  and  increase  subsidence  over 
the  surface  anticyclone  leading  to  stronger  surges.  Weak 
monsoon  surges  were  associated  with  near  zonal  flow  at 
higher  altitudes,  this  describes  the  upper-level  flow  used 
in  these  experiments.  The  inclusion  of  one  or  both  of  these 
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APPENDIX  A 


MODEL  PRIMITIVE  EQUATIONS 


Equation  of  State: 

po  “  RT 


Hydrostatic  Equation: 


11  „  _ 

3o 


no 


Horizontal  Momentum  Equations: 


n/dv\  +  fk  HV  +  if  (M)  -  ftl  5_(4>a)  <* 

\dt/  a  6a 


Continuity  Equation: 


«n  +  v  • <nv)  +  s_(-ia) 
St  6a 
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1st  Law  of  Thermodynamics: 
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LIST  OF  SYMBOLS 

radius  of  the  earth 

amplitude  of  the  perturbation 

specific  heat  at  constant  pressure 

ccriclis  parameter 

value  of  coriolis  parameter  at  45N 

force 

gravitational  acceleration 
wave  number 

pressure  (P  =  P  ♦  ?•  ) 

surface  pressure  at  45N,  1013.25  mb 

maximum  model  pressure  height,  125  mb 

gas  constant  for  air 

temperature 

surface  temperature  at  45N,  266K 

zonal  velocity  component  (0=0  +  O') 

maximum  velocity,  surface:  5  m/s,  upper:  65  m/s 

meridianal  velocity  component  (V  =  v  +  7') 

terrain  height 

specific  volume 

halfwidth  of  jet 

dry  adiabatic  lapse  rate,  6.5  X/lcm 

potential  temperature 

longitude 

3.14159 

normalized  pressure 
vertical  coordinate 

vertical  velocity  in  sigma  cooriiinate  system 
model  time 

model  time  at  which  terrain  growth  begins 
model  time  at  which  terrain  reaches  final  heigh 
latitude 
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<*>,  latitude  of  jet 

%  latitude  at  which  mountain  range  is  center 

♦  geopoter.ti  al,  also  current  terrain  height 

%  maximum  terrain  height,  3 000  m 

ft  angular  velocity  of  the  earth 
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